Spectrum of the U(l) staggered Dirac operator in four dimensions 
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We compare the low-lying spectrum of the staggered Dirac operator in the confining phase of 
compact U(l) gauge theory on the lattice to predictions of chiral random matrix theory. The 
small eigenvalues contribute to the chiral condensate similar as for the SU(2) and SU(3) gauge 
groups. Agreement with the chiral unitary ensemble is observed below the Thouless energy, which is 
extracted from the data and found to scale with the lattice size according to theoretical predictions. 
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I. INTRODUCTION 

In recent years the spectrum of the Dirac operator 
in QCD and related theories has been studied in great 
detail, in particular with regard to its relation to chi- 
ral random matrix theory (RMT) R] and, more recently, 
partially quenched chiral perturbation theory g. Both 
the distribution of the small eigenvalues and the spectral 
correlations in the bulk of the spectrum are described 
by universal functions that can be computed analytically 
in these theories. "Universal" in this context means in- 
dependent of dynamical details and only dependent on 
certain global symmetries (and their spontaneous break- 
ing). The spectral correlations are only universal below 
a certain limiting energy, which is called the Thouless 
energy |^ because of analog situations first studied for 
disordered mesoscopic systems. This picture has been 
verified numerically in great detail by lattice calculations 
for the gauge groups SU(2) and SU(3) in four and three 
dimensions and for the Schwinger model in two dimen- 
sions, see Ref. Q for a summary. 

In this paper we study the staggered lattice Dirac oper- 
ator in quenched U(l) gauge theory in four Euclidean di- 
mensions. The bulk spectral correlations of this operator 
have been investigated earlier in Ref. Q. Here, we con- 
centrate on the low-lying eigenvalues. We compare their 
distribution to predictions of chiral RMT and estimate 
the Thouless energy. Our study is not done merely for 
the sake of completeness. We are also motivated by the 
fact that, because of the different topological structure of 
U(l), the physics governing the small Dirac eigenvalues 
may be different from the non-Abelian case. 

The standard lattice action describing compact U(l) 
gauge theory in 4d is given by 



where /3 = I/s^j ^i ^ Ux,^ — e'xjp{i9x,fx) and 9p = 
Ox,ti + Qx+fi,v-Ox+i>,t,-Ox,i> for V T^ ^i. For p < (3^ ~ 1-01, 
the theory is in the confinement phase, exhibiting a mass 
gap and monopole excitations Q. For /3 > /3c, the the- 
ory is in the Coulomb phase with a massless photon m. 
There are many interesting questions concerning the or- 
der of the transition between the two phases and the 
possibility of a nontrivial continuum limit for (3 — > /3~ 
S . Only the confinement phase exhibits chiral symmetry 
breaking, which has been addressed in a number of recent 
numerical studies pHuj- In the strong-coupling limit 
/? — > 0, chiral symmetry breaking follows rigorously from 
infrared bounds |lj] and has also been calculated explic- 
itly 111]. The broken phase is characterized by a chiral 
condensate that is determined by the small eigenvalues 
of the Dirac operator according to the Banks-Casher re- 
lation |l|l. 

In Sec. O, we compute the Dirac spectrum in both 
phases and investigate in more detail the properties of 
the small Dirac eigenvalues in the confinement phase. 
Section III discusses the Thouless energy that limits the 
universal regime described by chiral RMT, and conclu- 
sions are drawn in Sec. tt\1. 



II. SMALL DIRAC EIGENVALUES 

The staggered Dirac operator is constructed from the 
gauge fields according to 



Dxy = ^Xl[^^(^)^'''^^)'^; 
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5{C/,}-/3^(l-cos0p) , 
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where a is the lattice spacing, which we shall set to unity 
in the following, and the 77^ are the staggered phases. For 
the purpose of comparing the spectrum of D to RMT 
predictions, we note that D is in the symmetry class of 



the chiral unitary ensemble (chUE) of RMT because it 
has complex matrix elements and no anti-unitary sym- 
metries. 

The nonzero eigenvalues of D come in pairs ±zA„ with 
Xn real. For convenience, we refer to the A„ as the eigen- 
values. The spectral density of the Dirac operator is given 
by 



pW = (T.^(^ 



An 



(3) 



where the average is over all gauge field configurations, 
weighted by exp(— 5*) in the quenched theory. If chiral 
symmetry is spontaneously broken, the vacuum is char- 
acterized by a nonzero order parameter, the chiral con- 
densate (V'^)- The Banks-Casher relation ||lj] states that 



E=|(V;V)| = lini lim TTp{e)/V , 



(4) 



where V is the four-volume. The order of the limits in 
this equation is important. Note that the condensate 
is due to an accumulation of Dirac eigenvalues close to 
A = 0. The Dirac operator can also have eigenvalues 
equal to zero, but this is not the case for the staggered 
Dirac operator at finite lattice spacing. 

If the Dirac spectrum corresponds to one of the RMT 
universality classes and supports a nonzero value of S, 
the distribution of the smallest Dirac eigenvalues is de- 
scribed by the microscopic spectral density |15 
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The quantity ps is a universal function that depends only 
on the number of massless (or very light) flavors N f and 
on the topological charge v^ which is equal to the num- 
ber of exact zero modes of D that are stable under small 
perturbations of the gauge field. The superscript {v) in 
Eq. (g) means that the average according to Eq. (||) is 
only over the configurations with topological charge equal 
\.o V. In our case, we have N j = since we study the 
quenched theory. Furthermore, we take v = {i because 
we are using staggered fermions which do not have ex- 
act zero modes at finite lattice spacing [llB]. This point 
has been discussed in Refs. [ p7|jlq ], and the only situa- 
tion where deviations from the result for v ~Q have been 
observed with staggered fermions is the Schwinger model 
in two dimensions at very weak coupling |19| . On the 
other hand, Neuberger's Overlap Dirac operator [gO[ al- 
lows for exact topological zero modes on the lattice, and 
lattice simulations with this operator indeed find agree- 
ment with the RMT predictions for t/ 7^ Hll. 

The microscopic spectral density can be computed an- 
alytically. The prediction of the chUE of RMT for this 
quantity is, for N f = i^ = 0, p^ 



Ps{z) = -{JI{Z)^JI{Z)\, 



(6) 
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FIG. 1. Spectral density of the staggered Dirac operator on 
identical scales in the confinement (left) and Coulomb (right) 
phase, respectively. The average is over 20 configurations. 
A nonzero value of E is supported only in the confinement 
phase. 



where J denotes the Bessel function. We also consider 
the distribution of the smallest eigenvalue of D for which 
the RMT result for iV/ = i^ = reads p) 



-P (Amin) 



(FS)2A 
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A comparison of lattice data with these predictions is 
only sensible if S > 0, i.e. if there is a sufficiently strong 
accumulation of small Dirac eigenvalues in the vicinity 
of A = 0. In Fig. H we have plotted the spectral density 
of the staggered Dirac operator in this region, computed 
on an 8"^ X 6 lattice for /3 = 0.9 (confinement phase) 
and f5 — \.\ (Coulomb phase), respectively. Clearly, a 
nonzero value of E is supported only in the confinement 
phase, and thus the following analysis will be done only 
for this phase. 

In non-Abelian gauge theories, the accumulation of the 
small Dirac eigenvalues is usually attributed to the pres- 
ence of instantons. The argument is that the degeneracy 
of the exactly zero eigenvalues in the field of isolated in- 
stantons is lifted by interactions, leading to eigenvalue 
repulsion and to a nonzero value of E. The topological 
structure of U(l) gauge theory in 4d is different, and ev- 
idence has been presented ||^ which suggests that mag- 
netic monopoles account for chiral symmetry breaking 
in the Abelian gauge theory. However, it is not quite 
clear whether the monopoles are really the driving mech- 
anism, or if disorder alone would be sufficient, because it 
is difficult to disentangle disorder and monopole effects 
convincingly. Neither the rigorous arguments P4 I nor 
the strong-coupling investigation ||l^ make use of any 
explicit mechanism. 

Let us turn to the analysis of our data. We have com- 
puted the eigenvalues of the Dirac operator on lattices of 
size 4**, 6^, and 8'^ x 6 using /3 = 0.9 in the confinement 
phase. To compare the data to the RMT predictions, we 
determine the parameter E in Eq. (0) by extrapolating 
the spectral density many level spacings away from zero 
to A = 0. This procedure is completely independent of 



V 



'x6 



config. 

10,000 

10,000 

3,745 



■^BC 



0.352(8) 
0.352(4) 
0.353(3) 



0.345(2) 
0.353(2) 
0.352(3) 



TABLE I. Summary of our simulations at /3 = 0.9. The 
parameter E was obtained by two different procedures as de- 
scribed in the text. 

RMT. As a check, we have also determined E via RMT: 
Using Eq. (^) , the expectation value of Amin is given by 

(A„,i„) = V^/{V^) , (8) 

which allows us to determine S from the numerical value 
of (Amin). Together with the numbers of configurations 
per parameter set, the values of S obtained from these 
two procedures are given in Table ||. The two values of S 
are in excellent agreement, except for the smallest lattice 
size, where the agreement is not perfect but still within 
error bars. 

In Fig. g, we have plotted the microscopic spectral den- 
sity and the distribution of the smallest eigenvalue for all 
three lattice sizes along with the predictions of the chUE 
of RMT. The lattice data for P(Ai„in) agree perfectly 
with Eq. (0). The microscopic spectral density (^) is also 
well described by RMT, but the agreement breaks down 
for large values of z, with a "critical" value of z that in- 
creases with lattice size. This is essentially the Thouless 
energy to which we now turn our attention. 



III. THOULESS ENERGY 

As mentioned earlier, the small Dirac eigenvalues are 
described by universal functions only for energies below 
a limiting scale, the Thouless energy. In QCD, this scale 
is determined by the requirement that the inverse mass 
of the pion is equal to the largest linear size of the box 
with volume V = L^^ x Lt, i.e., XjraT^ « max(Ls, Lt) 
||. This can be translated to E^ - fl/{Y.Ll) ||,|6|, 
where /^r is the pion decay constant and in our case we 
have Lg > Lt- A dimensionless estimate of the Thouless 
energy is obtained by expressing Ec in units of the mean 
level spacing at A = 0, given by A = l/p(0) ~ tt/{VJ^). 
This yields 

uc^^- -flL.U ■ (9) 

To test this prediction, we follow the lines of Ref. [M and 
construct the disconnected scalar susceptibility, defined 
on the lattice by 



Alatt 
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^ (iAfc -I- m){iXi -\- m) 



N 



■N\]^^^X, + m^ 



(10) 




FIG. 2. Microscopic spectral density (left) and distribution 
of the smallest eigenvalue (right) of the Dirac operator for 
three different lattice sizes. The histograms represent lattice 
data, and the solid lines are the RMT predictions. 



where m is a valence quark mass. The correspon dins 
RMT result for the quenched chUE with ly — reads |2q 



Xr'mt (") 



u'[I^{u)-lt{u)][Kt{u)-K^{u)], (11) 



where u — mVY, and / and K are modified Bessel and 
Hankel functions. The quantity xf^tt should be described 
by Eq. (pT| ) for u < Uc- The dimensionless Thouless 
energy can be extracted by inspecting the ratio |E7| 



ratio =(xfi\^-X^^%)/x 



RMT 



(12) 



This quantity should be zero for u < Uc and deviate 
from zero for u > u^.. The data for this ratio computed 
at /? = 0.9 for our three lattice sizes are shown in Fig. ||. 
Consider first the left plot. It is clear that Uc increases 
with increasing lattice size. To test the scaling predicted 
by Eq. (Bh, the same data are shown in the right plot, 
but now plotted versus u/{LsLt). The data for different 
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FIG. 3. The ratio of Eq. (|]J) plotted versus u (left) and 
u/{LsLt) (right). The deviations of the ratio from zero for 
very small values of u are well-understood artifacts of the 
finite lattice size and finite statistics [|27[. 



lattice sizes now fall on the same curve, confirming the 
predicted scaling behavior of the Thouless energy. 



IV. DISCUSSION 

We have shown that in the confinement phase of com- 
pact U(l) gauge theory on the lattice, the distribution 
of the small Dirac eigenvalues is described by universal 
functions that can be computed in chiral RMT. The lim- 
iting energy above which non-universal behavior emerges 
scales with the lattice size as expected. 

The origin of the small eigenvalues in U(l) gauge the- 
ory deserves further attention. The question may be less 
about a mechanism in the strong-coupling limit, where 
it appears that the disorder of the gauge fields could be 
sufficient. The important question is about a mecha- 
nism that could sustain chiral symmetry breaking for 
a large correlation length, eventually leading to a con- 
fined QED continuum theory for (3 -^ (3~ . There, U(l) 
monopoles could play a crucial role. Instantons appear 
to provide such a mechanism for SU(2) and SU(3) non- 
Abelian gauge theories. 

In the Coulomb phase, chiral symmetry is restored, so 
the "critical" value (3'^ for the chiral phase transition can- 
not be larger than [3c- However, we know of no strict ar- 
gument that confinement implies chiral symmetry break- 
ing, so it is possible, at least in principle, that (3'^ < Pc- 
(In supersymmetric theories, one can have confinement 
without chiral symmetry breaking |Q.) Because the chi- 
ral condensate is directly related to the distribution of 
the small eigenvalues, the chiral phase transition can be 
studied by observing the distribution of, say, the small- 
est positive eigenvalue for /3 — > /3~ Ig^l- This is another 
reason why it would be interesting to study this limit in 
future work. 
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